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ABSTRACT 

Recently, Hennebelle and Chabrier (2008) derived an analytical theory for 
the mass spectrum of non self-gravitating clumps associated with overdensities 
in molecular clouds and for the initial mass function of gravitationally bound 
prestellar cores, as produced by the turbulent collapse of the cloud. In this 
companion paper, we examine the effects of the non-isothermality of the flow, of 
the turbulence forcing and of local fluctuation of the velocity dispersion, on the 
mass function. In particular, we investigate the influence of a polytropic equation 
of state and of the effective adiabatic exponent 7 and find that it has a drastic 
influence on the low mass part of the IMF. We also consider a barotropic equation 
of state (i.e. a piecewise polytropic eos) that mimics the thermal behaviour of 
the molecular gas and compare the prediction of our theory with the results 
of numerical simulations and with the observationally-derived IMF, for cloud 
parameters which satisfy Larson's type relations. We find that for clouds whose 
density is, at all scales, almost an order of magnitude larger than the density 
inferred for the CO clumps in the Galaxy, a good agreement is obtained between 
the theory and the observed IMF, suggesting that star formation preferentially 
occurs in high density environments. We derive an analytical expression for the 
IMF which generalizes the expression previously obtained for the isothermal case. 
This easy-to-implement analytical IMF should serve as a template to compare 
observational or numerical results with the theory. 

Subject headings: stars: formation — stars: mass function — ISM: clouds - 
physical processes: turbulence 
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1. Introduction 

Understanding star formation is one of the greatest challenges of astrophysics. Amongst 
many outstanding questions, finding out what determines the stellar initial mass function 
(IMF) first observationally derived by Salpeter (1955) and later on by many others (e.g. 
Kroupa 2002, Chabrier 2003ab, 2005), constitutes a key issue. Various observational works 
have recently suggested that the IMF could be inherited from the conditions prevailing in 
the molecular clouds. In particular the core mass function (CMF) has been found to have 
a shape very similar to the IMF (Motte et al. 1998, Testi & Sargent 1998, Johnstone et al. 
2000, Alves et al. 2007, Nutter & Ward-Thompson 2007, Simpson et al. 2008, Enoch et 
al. 2008), although shifted towards larger masses by a factor of about 2-3. Although this 
behaviour remains to be confirmed unambiguously, it suggests that the IMF is rooted in the 
physics of the turbulent self-gravitating gas which very likely determines the core formation. 
This interpretation has been questioned by Clark et al. (2007) who point out that since the 
lifetime of the more massive cores is likely to be longer than the lifetime of the less massive 
ones, the core mass function should be indeed flatter in order to lead to the IMF. This issue 
is investigated in appendix C where a simple solution is proposed. The relationship between 
the CMF and the IMF has also been recently discussed by Goodwin et al. (2008) and 
Smith et al. (2008). Theoretically, IMFs very similar to the one inferred from observations 
have been obtained from hydrodynamical numerical simulations of supersonic, isothermal 
self-gravitating gas (e.g. Klessen 2001, Bate & Bonnell 2005). An analytical approach of 
the IMF based on such a turbulence-driven star formation picture, although ignoring the 
turbulent support, has been proposed by Padoan et al. (1997) and Padoan & Nordlund 
(2002). These authors conclude that in the absence of magnetic field, a too stiff IMF is 
obtained and that one needs an at least weakly magnetized gas to obtain a CMF/IMF with 
a slope close to the Salpeter's value, contradicting the results of Bate & Bonnell (2005) and 
Jappsen et al. (2005). 

Recently Hennebelle & Chabrier (2008, hereafter paper I), have extended the formalism 
derived in cosmology by Press & Schechter (1974) to predict the galaxy formation mass 
spectrum, to the context of star formation, and have derived a new analytical theory of the 
IMF. These authors show that the IMF includes two contributions, namely a power-law at 
large scales, and a lognormal contribution which yields a truncation at both large and small 
scales. The first contribution stems from the scale introduced for gravitational instability, the 
Jeans scale in the present context, while the second contribution arises from the statistical 
counting of overdense structures exceeding some density threshold, similar to the Press- 
Schechter cosmological case. This theory is based on the same gravo-turbulent picture as 
the one proposed by Padoan et al. (1997) but, among other differences, explicitly includes 
the turbulent support, emphasizing the dual role played by turbulence, which leads to the 
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enhancement of (collapsing) overdensities on one hand and to a non-thermal support of the 
cores on the other hand. The theory predicts a CMF/IMF very close to the observational 
one, even in the hydrodynamical case, without the need to invoke magnetic field, and relates 
explicitly the slope of the high-mass power-law part of the IMF to the power spectrum 
index of turbulence. Within the same formalism, when invoking only the density threshold 
condition, the mass spectrum of non self-gravitating CO clumps, identified as overdensities in 
molecular clouds, is also derived and is found to be in good agreement with the mass spectrum 
inferred observationally (e.g. Heithausen et al. 1998), i.e. dN/dM = J\f(M) oc M~ x with 
x ~ 1.6- 1.8. 

In the present paper, we extend various aspects of this theory, exploring the depen- 
dence of the IMF upon various thermal and dynamical properties of the flow, namely non- 
isothermality, turbulence driving and local fluctuations of the velocity dispersion. As will 
be seen, both aspects are playing a significant quantitative role although not modifying 
qualitatively the main features of the theory. 

The paper is organized as follows: in §2, we summarize our theory and in §3 we inves- 
tigate the influence of the velocity statistics on this theory. In §4 we consider a polytropic 
equation of state while in §5 we study the influence of a barotropic equation of state that 
mimics the temperature distribution of the molecular gas. In §6 we compare the theory 
prediction with the results of numerical simulations and with the IMF inferred by Chabrier 
(2003, 2005), determining the set of cloud parameters for which good agreement is obtained. 
Section 7 concludes the paper. 

2. Summary of the theory 

The theory is based on an extension of the statistical Press & Schechter (1974, PS) 
formalism (see also Bond et al. 1991, Padmanabhan 1993, Jedamzik 1995, Audit et al. 
1997), developed in cosmology, to the context of star formation, characterized by a non- 
uniform (lognormal) underlying density field. The principle of the method is the following. 
First, the density field is smoothed at scale, R, using a window function (which in our case 
is the sharply truncated in k-space window). Then, the mass contained in areas which at 
scale, R, satisfy a specified density criterion is counted. This mass is equal to the mass that 
will end up in structures of mass larger than a scale dependent critical mass. Taking the 
derivative with respect to R of the corresponding equation, we obtain the mass spectrum. 

In our theory, non self-gravitating CO clumps are identified as local overdensities, 5 = 
log(p/ p), exceeding a simple, constant density threshold 5 C , while in the case of gravitationally- 
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bound prestellar cores, the density threshold is given by the requirement that a fluctuation 
contains one local (thermal or turbulent) Jeans mass. Whereas pure thermal support yields 
a CMF at large scales steeper than the Salpeter one, dN/dlog M oc M~ 135 , turbulent dis- 
persion leads to the correct slope, for the index of turbulence found in observations or in 
numerical simulations of isothermal compressible turbulence (e.g. Kritsuk et al. 2007). Tur- 
bulent support (e.g. Bonnazola et al. 1987, Vazquez- Semadeni & Gazol 1995) enters the 
theory through an effective Mach number: 



defined as the ratio of non-thermal vs sound speed at the mean Jeans scale A° (and not at 
the local Jeans length), whereas the usual Mach number, Ai, represents the same quantity 
at the scale of the turbulence injection scale, L iy assumed to be the characteristic size of the 
system: 

M = (2) 

Since the global Mach number, Ai, broadens the density PDF, it describes the trends of 
supersonic turbulence to promote star formation by creating new collapse seeds while on 
the contrary, the effect described by Al* is to prevent the gravitational collapse through the 
turbulent dispersion. Here, C s = (kT '/ ' ^m H ) 1 ' 2 « 0.2 (/i/2.0)~ 1/2 (T/10K) 1 / 2 kms" 1 denotes 
the thermal sound speed, ji is the mean molecular weight, and V is the (scale-dependent) 
rms velocity which obeys the Larson relation (Larson 1981): 

(Ui) 2 ) = ^x(^) 2 ', (3) 

with Vo — 1 kms -1 and rj ~ 0.4-0.5. The various symbols and notations used in the paper 
are defined in Table 1 of paper I. 

The scale-dependence for the velocity dispersion associated with supersonic turbulence 
adopted in our theory is given by (see paper I for details): 
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with 



^ = ln(l + 6 2 .M 2 ) 
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where 5 is the Fourier transform of 5 = log(p/p) and b is a constant obtained from numerical 
simulations^. In this expression, the Mach number A4 can be either a hydrodynamical or 
an Alfvenic Mach number (Vazquez-Semadeni 1994, Padoan et al. 1997, Passot & Vazquez- 
Semadeni 1998, Ostriker et al. 2001, Kritsuk et al. 2007). As mentioned in paper I, 
in case of supersonic turbulence, the spectral index, n', of log(p) calculated in isothermal 
hydrodynamical and MHD simulations is found to be close to the n = 11/3 value obtained 
in incompressible turbulence for the velocity field (Beresnyak et al. 2005, Federrath et al. 
2008). In our approach, the star- forming clumps issued from these large-scale turbulent 
motions are identified with over-densities 5 = log(p/p) and the mass associated with these 
over-densities at scale R is 

M ~C m R 3 p = C m R 3 pe 5 , (6) 

where C m is a geometrical coefficient of the order unity which depends on the window 
function, as discussed by Lacey & Cole (1994). For a top hat function, it may be simply 
equal to 4n/3 whereas for the sharply truncated function in k-space, these authors suggest 
the value 6% 2 . This gives for the smoothing length, R = [Mj (C m p)) 1//3 exp(— 5/3). Unlike in 
the PS formalism, the mass depends not only on the spatial scale, R, but also on the variable 
5, a consequence of the non-uniform, lognormal density fluctuation distribution produced by 
turbulence. 

The mass spectrum of the number-density, Af(M) = dN/dM, of non self-gravitating 
clumps defined by a density threshold 5 C is found to be (see also Hennebelle & Audit 2007): 



where M = C m pLf is the whole mass contained within a volume C m L\ and p c = pe Sc . 

Within the same theoretical formalism but with a density threshold which now depends 
on density and thus on the scale, namely the Jeans stability criterion, the IMF of self- 
gravitating cores is found to be: 

x \n(M/B?) 

Af(M) = 2K i 1 + < 1 - " ) - M - fl2 '' x ( ™ ) " x -P<- ff2 /8) 



R e [1 + (2rj + l)MlR 2r i] \R 3 J V2ira 



^^Note that in paper I, we wrote b instead of b 2 in eqn. ([S]), but the quantitative value was consistent with 
the widely used value b 2 « 0.25 (e.g. Padoan et al. 1997) 
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where R = R/X°j, M = M/M% = R(1+ M 2 R 2r >), 6 C R = In {(1 + MlR 2l1 )/R 2 Y Af = p/M°j 
and Mj, X°j denote the usual thermal Jeans mass and pseudo Jeans length|§, respectively: 

M ° = aj -j2= « 1.0 aj (JL)3/2(_^_)-i/2 ( n J -V2 M (9) 
y/CPp V 10K ; V 2.33 ; l 10 4 cm-3 ; ° V; 

1 /3 

A o = f^L) -^^o.la 1 /^— ) 1 ^(JL)-i/2( " J " 1 / 2 pc (10) 
J \C„J v 7 ^ 10K 2 - 33 V 10 4 cm- 3; P v ; 

where aj is a dimensionless geometrical factor of order unity. Taking for example the stan- 
dard definition of the Jeans mass, as the mass enclosed in a sphere of diameter equal to the 
Jeans length, we get aj = 7r 5 / 2 /6 and in this case the pseudo Jeans length is just equal to 
the Jeans length divided by 2. Note, however, that the exact derivation of aj is not very 
accurately determined since it is related to the criterion for a mass of gas to collapse. In 
particular, since the density distribution is not homogeneous and since turbulent dispersion 
or support also enter the calculations, no precise value can be inferred. As seen in § 16. 1\ 
the above mentioned estimate leads to good agreement with the simulations, even though a 
more refined calibration might be required at some stage. 

The exponent rj in eqns. (pQ) and (j3J) is related to the aforementioned (3D) index of 
turbulence n by 

n — 3 . . 

■ ( U ) 

As mentioned above and in paper I, strictly speaking, the index of turbulence n which appears 
in this expression, which is related to the velocity power spectrum, is not necessarily the same 
as the one which enters eqn. (j3j), nf, which is related to the power spectrum of log p. However, 
as mentioned above, numerical simulations seem to find that both indexes are similar. Recent 
high resolution simulations of non-magnetized isothermal supersonic turbulence (Kritsuk et 
al. 2007) yield n ~ 3.8-3.9, i.e. r] ~ 0.4-0.45. 

We recall (see paper I) that the present theory yields an analytical relationship be- 
tween the high-mass power-law index of the aforementioned clump and core mass functions, 
dN jdM oc M~( 1+x \ and the turbulence power spectrum indexes, such as: 

n' 

x = 2 for the non self — gravitating clumps (12) 

3 

n + 1 

x ~ for the self — gravitating cores, (13) 

2n-4 6 & y 1 



2 Strictly speaking the thermal Jeans length is y/nC s /y/Up, the pseudo Jeans length differs from it by a 
constant factor. 
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which, for the second relation, yields the Salpeter value, 1 + x ~ 2.33 for n ~ 3.8, although, 
strictly speaking, the exact relation involves a second correcting term (see §6.2.2 below). 

Note that, as evident from eq. (jlj) of the present paper and eq. (33) of paper I (which 
entails the term da/dR arising from the calculation of dP^/dR), the value n' = 3 appears to 
be critical. If n' < 3, small scale density fluctuations are dominant while the contrary is true 
when n' > 3. As recalled above, tridimensional turbulent flows produce powerspectra whose 
index at scale smaller than the injection scale, is typically of the order of 11/3 and therefore 
larger than 3. However, driving the turbulence at small scales, would produce powerspectra 
with different scale dependence and would therefore have serious consequences on the mass 
spectra since it would modify the dependence of o on R. This is consistent with the results 
of Klessen (2001) regarding the strong influence of the driving scale on the mass spectrum. 

3. Influence of the statistics of the turbulence 

3.1. Turbulence forcing 

The value b ~ 0.5 in eqn.(15]) has been derived from numerical simulations and seems to 
apply to both hydrodynamical and MHD simulations (Padoan et al. 1997, Li et al. 2008). 
However, these simulations include only solenoidal (divergence-free) modes of turbulence 
forcing. Simulations including the contribution of the compressive (curl-free) modes in the 
forcing (Federrath et al. 2008, Schmidt et al. 2009) lead to a value of about b ~ 1 for 
the purely compressible forcing case, b ~ 1/3 for the 3D case with only solenoidal forcing, 
and values between ~ 1/3 and 1 when both types of modes contribute, as expected for 
the turbulence taking place in the interstellar medium. It should be noted, however, that 
Federrath et al. infer PDFs that are non perfectly Gaussian, particularly at high density. 
Investigating the influence of such a distribution remains an open issue. Figure [T] portrays 
the theoretical CMF/IMF (eqn.flBJ) for b = 0.5, b = 0.75 and 6 = 1. While the first value 
corresponds to about equipartition of solenoidal and compressive modes, 2:1 in 3D (equation 
5 of Federrath et al. 2008 with £ = 1/3), the intermediate one corresponds roughly to 
equipartition between total solenoidal and compressive energy, i.e. the compressive modes 
have twice as much energy as the solenoidal modes. Not surprisingly, including compressive 
forcing yields a broader PDF for turbulence and thus increases the number of small-scale 
overdense regions. This shows that, as mentioned by Federrath et al. (2008), the statistics 
(PDF) generated by compressible turbulence depends not only on the Mach number but 
also on the forcing and on the relative contributions of rotationally and compressibly driven 
modes. This double dependence enters through the values of M. and b in eqn.([5]) in the 
present theory. Indeed, since a 2 = ln(l + b 2 Ai 2 ), changing b is equivalent to changing M. 
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in bAi. Thus, b = 1 is equivalent to the b = 0.5 case but for a Mach number twice smaller. 
Since, as mentioned in paper I (fig. 5), good agreement between the present theory with 
b = 0.5 and the Chabrier's IMF requires Mach numbers of about 12 (mostly for the low mass 
end of the IMF), this implies that a Mach number of ~ 6 will be sufficient if b ~ 1 to get 
a good agreement. Except if otherwise specified, the calculations in the following sections 
have been conducted with a value b = 0.5. 

To end this section, we want to stress the following issue. Remembering that the 
parameter M.* represents the turbulent support at the Jeans length, it is likely that, when 
considering flows dominated by compressible motions, its value should be reduced, since 
in that case a significant fraction of the modes are not quenching the collapse but, on the 
opposite, may tend to promote it. 



As recalled in §[2J the present theory crucially depends on the parameter .M* (eqn. (JT])) 
which represents the velocity dispersion at the Jeans length characteristic of the cloud's mean 
density. As seen from eqn. (JBJ), a single averaged value of .M* has been used so far. In reality, 
one expects the velocity dispersion to undergo fluctuations within clouds. The dispersion of 
the local Mach number has been computed in numerical simulations (e.g. Kritsuk et al. 2007, 
fig. 4). A somehow broad distribution has been obtained showing no obvious correlation with 
the gas density. The shape of this distribution appears to be complex though clearly peaked 
around a mean value. For sake of simplicity, we will assume in this work that the distribution 
of the rms velocity, P rms , is lognormal with a width that we estimate from the numerical 
simulations. We write, with X = ln(.M*): 



where lnA^* denotes the mean value of ln.M*. From fig. 7 of Kritsuk et al. (2007), we 
estimate that a rms ~ 0.3 — 0.5, although the distribution is not exactly lognormal. Given 
the lack of knowledge of P rmsi we further assume that a rms does not depend on the scale, 
i.e. that the dispersion of the rms velocity is the same at all scales. 

In order to obtain the mean core mass spectrum, J\f(M, In A4*), we must now integrate 
the core mass spectrum obtained for a single value of .M*, Af(M, In .M*) (eqns. (33) and Bl 



3.2. Fluctuation of the rms velocity within the cloud 



3.2.1. Principles 




(14) 
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of paper I and eqn. (jSJ) of this paper), over the A4* distribution: 

POO 

AT(M,WW,) = / P rTM (X,KM m )JJ'(M,X)dX (15) 



3.2.2. Results 

Figure [2] shows the core distribution for various values of a rms and for M. = 6 and 
ln.M* = In 2/2. Taking into account a distribution of M.* has a clear though modest effect 
on the mass spectrum by modifying the resulting core distribution in various ways. While 
the amount of small-mass objects is almost unchanged, more stars form in the intermediate 
mass range whereas less massive stars are produced, resulting in a slightly stiffer CMF/IMF 
with increasing a rms . This is reminiscent of the trend displayed in fig. 1 of paper I which 
shows that small values of .M* tend to produce more intermediate mass objects and less 
massive stars. Since the core mass function obtained by taking into account the dispersion 
of Ai* resembles the one obtained with a unique but lower value of .M*, this implies that the 
low part of the .M* distribution has more influence than the high part of the distribution. 
However, since the overall effect of introducing a distribution of M.* values rather than a fix 
value appears to be modest, a constant .M* will be considered in the rest of the calculations. 



4. Polytropic equation of state 

It is well established that including a detailed treatment of the thermal properties of the 
gas in the molecular cloud leads to an equation of state, P = Kp 1 ', softer than isothermal, 
with 7 « 0.7-0.8 (Larson 1985, Scalo et al. 1998, Glover & Mac Low 2007). In this section, 
we explore the consequences of the departure of the gas from isothermality on the CMF/IMF. 



4.1. Thermal criterion with a polytropic equation of state 

We first investigate the case of a simple polytropic equation of state (eos). The eos 
reads, assuming the flow behaves perfect gas: 

P = K p< = — p T(p), (16) 
pm H 

which leads for the constant K: 

K( 1 ) = K = —p 1 -iT(p). (17) 
pm H 
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log(M/M„) 

Fig. I. — Core mass function for a Mach number M = 6 and Ml = 2. Various values 
of the coefficient b are shown, namely b=0.5, 0.75 and 1 from right to left. The solid and 
dotted lines in the upper right corner illustrate the behaviour of a Salpeter distribution, 
dN '/dlogM oc M~ L35 and a slope dM/ dlogM oc M~ L0 , respectively. 




-3-2-10 1 2 

log(M/M ) 

Fig. 2. — Core mass spectrum for various values of <J rms , the width of the .M* dispersion. 
While the dotted line shows the case a rms = 0, the solid lines display the cases a rms = 0, 
0.2, 0.4 and 0.8. 
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Conversely, the temperature now depends on the density as: 

T = To(^) 1 " 7 , (18) 
P 

where T and p = p denote the average initial temperature and density of the cloud. The 
thermal sound speed now depends on the local density (e.g. McKee et al. 1993, Vazquez- 
Semadeni et al. 1996) as well as: 



cj~) = J^ = ^(7)(^)^=c<«( 7 K-- no) 

withC?( 7 ) = M=( 1 K)^- p ^ = (^LT )^p 1 -^, (20) 
V op P™h 

where 5 = ln(p/p) is the density threshold for overdense structures. 

Accordingly, the mean thermal Jeans mass and pseudo Jeans length now become: 

M°'( 7 ) = aj^^l = M«^\ (21) 
V G P 



where Mj and Aj are given by eqns.(!9|)- (|T0|) . The local thermal Jeans mass and pseudo Jeans 
length thus read: 

Mj = a J ^ = M°'( 7 )(^)^- 2 = M°'( 7 )e^- 2 ) 5 , (22) 
V G P p 

A j = (^)V3 = A y (7)e (ir-i}«. 

The conditions for thermal collapse (eqns. (19) of paper I) now become: 

M > Mft = M°jR { ^\ (23) 

5>5 C R = -^ln(£), 
7-2 

where R = R/X°j. 
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4.2. General criteria 

Proceeding as in paper I, i.e. denning an effective sound speed and Jeans mass, the 
threshold condition for gravitational collapse in the general case where both thermal and 
non-thermal support contribute (eqns. (28)-(29) of paper I) now reads: 

M>Mj = aj [ s{1) + 3^ mgi_. (24) 



With R = R/X°J(j) and M = M/Mj"( 7 ), this leads to 

M>M C R = R\(~Y +MlR 2v 




(25) 

where Ml = MHj) = § {V /C°(j)) 2 x (\°J/lpc) 2r >. Unfortunately, eqn. ([25]) cannot be 
analytically solved and M R must be determined numerically: 

W R = +MIrA. (26) 

Although resolving eqn. ( 1261) is not difficult, it is convenient to have an analytic expression 
such as eqn. (JSJ) (i.e. eqn. (44) of paper I), for the general 7 ^ 1 case. In appendix A, we 
we derive an accurate and analytic approximation of eqn. (|26|) for such a general IMF. This 
analytic expression can serve as a quick- and- easy first approach to compare IMFs derived 
from either observations or numerical simulations with the present theory. 



4.3. Core mass function with a polytropic equation of state 

Before we derive the core mass function, it is necessary to discuss the density PDF. As 
recalled in §2, the distribution of the logarithmic density fluctuations produced by isothermal 
turbulence is a Gaussian with a width given by eqns. (j4]) and (jSJ). It seems obvious from these 
equations that a more compressible flow (7 < 1) will increase the variance of the density 
PDF for a given injection scale and a given velocity dispersion. 

Indeed, the relation stated by eqn. (jSJ) is valid only in the isothermal case and no 
extension has been proposed when 7 7^ 1 in 3D. The density PDF itself has not been studied 
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and could be different from a lognormal distribution, as discussed in Passot & Vazquez- 
Semadeni (1998) in the case of ID flows and, though with a limited numerical resolution, by 
Li et al. (2003) in 3D. We also recall that Federrath et al. (2008) find substantial departures 
from a lognormal distribution when forcing in the compressible modes, even in the isothermal 
case, with a significantly broader PDF. Since this issue appears to be unsolved yet, in this 
work we will simply use eqn. (j3j) even when 7^1. By doing so, we probably, for a given 
Mach number, underestimate the width of the density PDF when 7 < 1 and overestimate 
it when 7 > 1. Since, in the following we use a fiducial Mach number M. = 6 in our 
calculations, this implies that, rather than a fixed Mach number, we assume a fixed PDF 
width, a, which may correspond to Mach numbers different from M. = 6 as 7 varies. Given 
the lack of knowledge on the nature of turbulence in non-isothermal flows, it seems difficult 
to go further at this stage. 

The mass spectrum, Af, is given by eqn. (33) of paper I, leading to (neglecting the 
large-scale second term in the complete derivation, see §5.1.2, 5.4 and appendix B of paper 

I): 

,r^r, P dR d8 c R 1 / (5 C R ) 2 5% a 2 \ , , 

Js[(M%) = ^ — ^ £ ; exp + , 27 

R MjM R dM R dR V^tto^ 2a 2 2 8 J' 



where 



^ = ln(S|, (28) 



dM R 



dR 1 - (7- l)(M R p~ 2 R^ 



((4 - 3 7 ) (M^) 7_1 i? 3 ~ 37 + (2 V + l)MlR 2 ^ , (29) 



dd _ dM c R 1 3 
dR dR M C R R 



(30) 



It is easily seen that for 7 = 1, the isothermal limit investigated in paper I is recovered, 
while the case of pure thermal support is recovered for M.+ = 0. 



4.4. Peak position for a polytropic equation of state 



To obtain the peak position of the CMF, we proceed as in paper I, i.e. we derive the 
core mass function, A/", of the purely thermal case (i.e. M.* = 0) with respect to the mass, 
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and calculate the mass given by the conditions dN jdM = 0. Taking into account thermal 
support only is a valid approximation since, as shown in fig. 1 of paper I, .M* has a weak 
influence on the peak position. 



M peak T , / (4-3 7 )(-9 + 6 7 )a 2 . , , ,, 

- M peak = exp ( ] (31) 



This yields: 



M — p«k y 4 

= (l + b 2 M 2 ) * 

This equation shows that the peak position strongly varies with 7 when 1 + b 2 M. 2 > 1. For 
example, for 1 + b 2 M 2 = 10, 7 = 1 leads to M pe ak = 10" 3/4 ~ 1/6 while 7 = 0.7 leads to 
Mp ea k — 1 /10 2 " 3 ~ 1/200. A small decrease of 7 thus leads to a large variation of the peak 
of the CMF. This stems for the stiff variation of the critical mass and density (eqns. (T20]) . 
(12"4"|) and (1251) ) with 7. As expected, when 7 — > 4/3, the peak of the CMF tends toward 
M peak = 1, i.e. M peak -> M». 

It is also interesting to compute the powerlaw index of the mass spectrum when turbu- 
lence is weak, (M* ~ 0). In this regime, we find that 

M{M) oc M"^ = M" (1+x) (32) 

This shows that the larger 7, the stiffer the core mass spectrum. For any reasonable value of 
7, the value of 1+x remains larger than l + (n+l)/(2n — 4) ~ 2.33 for n ~ 3.8 (see eqn. (j43l) 
and paper I), a consequence of the presence of turbulent support in the latter expression. 
Only when 7 < 0.2, does 1 + x become smaller than the Salpeter's value of ~ 2.3. 



4.5. Results for a polytropic equation of state 

Figure [3] compares the CMF/IMF obtained for various values of 7, namely 7 =1.3, 
1.1, 1, 0.9, 0.8, 0.7 for a typical Mach number M. — 6 and Ml = 2%. Varying 7 has a 
drastic influence on the core mass function. As expected from eqn. (131]) . more compressible 
turbulent flows (7 < 1) produce more small-scale overdense collapsing structures for a given 
a than the isothermal gas. For 7 < 0.8, the number of small objects increases dramatically. 
On the contrary, the number of small mass objects decreases abruptly when 7 > 1. We also 
see that the flow will not produce collapsing structures for 7> 1.3. These results arise from 



3 We recall that we assume a constant width a of the lognormal distribution, which would correspond to 
M. = 6 in the isothermal case. 
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the density dependence of the local Jeans mass, Mj oc p2^~2.) (see eqn. (1221) ) . which implies 
that the Jeans mass increases with density for 7 > 4/3, leading the collapse to choke itself, 
as pointed out in hydrodynamical simulations by Larson (1985), Vazquez-Semadeni et al. 
(1996) and Li et al. (2003). 

On the other hand, with the chosen value }A\ = 2, the number of large mass objects is 
almost unchanged when 7 varies. This stems from the fact that, as long as 7 > 0.2, turbulent 
support is dominant for large masses, as expected from eqn. (I3"2l . 



5. Barotropic equation of state 

As shown in the previous section, changing 7 has a drastic influence on the mass function. 
In particular, small values of 7 tend to produce more small mass objects which form at high 
density in the very cold gas. However, the thermal behaviour of the gas in molecular clouds is 
not expected to be well described by a single powerlaw (e.g. Larson 1985). Various authors 
(e.g. Glover & Mac Low 2007) found that for densities smaller than n ~ 10 5 cm -3 , i.e. 
p ~ 10~ 18 gem -3 , the effective polytropic exponent is about 7 ~ 0.7, whereas at higher 
densities 7 ~ 1 (Masunaga & Inutsuka 2000). As seen from fig. [31 the CMF obtained with 
7 = 0.7 is rather different from the one obtained for 7 = 1. Thus, taking into account the 
variations of 7 appears to be important. 



5.1. Core mass function with a barotropic equation of state 

Jappsen et al. (2005) adopt the following prescription for the temperature variation in 
the clouds: 



T = a l p'^ 1 p<p 



^crit 



crit 



T = a 2 p^~ l p>p 

where the critical density, p crit ~ 10~ 18 gem -3 corresponds to n crit 
71 =0.7 and 72 = 1.1. In our work, we find it more convenient to write 



(33) 

2.5 x 10 5 cm" 3 , 



{c o )2( P )72 -i 



m-i l/m 



\2 



(71— l)m 



+ 



r^O \ 2m / \ (72-l)m 
°s,2 \ IP 



l/m 



(34) 
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Fig. 3. — Core mass function for M. = 6 and Ail = 2 and various values of the polytropic 
exponent 7, namely 7 =1.3, 1.1, 1, 0.9, 0.8 and 0.7. 




Fig. 4. — Temperature distribution as stated by eqn. (15111 for m equals to m — 1, 2 and 3 
(from top to bottom). The dotted and dashed lines correspond to the prescription of Jappsen 
et al. (2005). 



10 20 30 40 50 



Fig. 5. — Peak position of the CMF as a function of 1 + b 2 Ai 2 . Solid lines display the peak 
position (eqn. for 7 2 = 1.1 and for various values of p crtt /p =1, 10, 10 2 and 10 3 

(from top to bottom). The dashed line corresponds to the peak position for ji — 0.7 
(eqn. EIJ) whereas the dotted line is for 71 = 1. 




log(M/M ) 

Fig. 6. — Core mass function for M. = 6 and M. 2 = 2, for various values of the constant ratio 
p CTlt /p, namely 10, 10 2 and 10 3 , from right to left. The solid lines correspond to the case 
m = 1 while the dashed lines correspond to m — 3. The dotted line displays the isothermal 
case. 
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where m is a real number of order unity and where C° = Cs(t) is given by eqn. fTSO]) . The 
exact dependence of the temperature upon density is not well known and may vary from 
place to place. Figure H] shows the temperature distribution as a function of the density. The 
two dash/dot-lines correspond to the prescription of Jappsen et al. (2005). In the following, 
we investigate how the CMF depends on the critical density p cnt and on the parameter m 
which, as seen from fig. HJ has some impact on the temperature variation. Low values of 
m correspond to a smooth transition between the two regimes whereas larger values lead to 
a more abrupt transition. As can be seen from fig. 2 of Larson (1985), our knowledge of 
the temperature distribution is not very accurate and it is unclear what would be the most 
appropriate choice for m. 

Since at the critical density, the two contributions are equal, we have the relation: 

(„crit \ 7i - 1 / „crit \ 72 - 1 

—) =( c y 2 (V) ■ (35) 



' C° \ 2 / n crit\ 71-72 

*- = ( & ) = { V ) ■ (36) 



leading to 



ciJ V p 

The value of K CT - lt depends on the cloud's mean density n and on the critical density, n cnt . 
Since molecular clouds have an average density of about n ~ 100 cm -3 , we find that the 
aforementioned value of n cnt leads to i^ C rit — 0.06. Taking a higher value, n ~ 10 3 cm -3 , 
which could be more representative of denser environments, yields K cvlt ~ 0.14. 

With these expressions, eqns. (|2"B"|) and (|2"9~|) now become: 

M C R = R (A 1 / 171 + MlR 2ri ) , (37) 



— \ (7i-l) m / — - \ (72-l)m 



and 



dM c R B 

~dR ~ C' _ _ 

with B = A^ + (3-3 7l )(^) m(7l - 1) A^- 1 + (ir crit ) m (3 - 3 72 )(^) m(72 - 1) A^ 

R 3 R 3 

+ (2rj + l)MlR 2 \ 
C = l-A^x 



(38) 



( 7l - l)^(S)^- 1 ) + (K CTit T^2 - ) . 

M C R R 3 M C R R 3 
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while eqns. ( I2"7|) . (1281) and fl30l) remain unchanged. 



5.2. Peak position for a barotropic equation of state 

Finding an exact expression for the position of the peak with a barotropic equation of 
state is not possible. However, one can estimate its value in two limiting cases. First, if 
the contribution of the second term to the thermal pressure (eqn. (|34"|) ) is negligible, i.e. 
K CT a <C 1, the eos is nearly polytropic and the peak position is given by eqn. fl3T|) . Second, 
if for the densities corresponding to the mass of the peak position, ~ M pea k/-R 3 , the second 
contribution of eqn. f[3~4"|) is dominant, then in the thermally dominated regime, we have with 
eqn. ([38]): 

. 1 ~ 4-372 

(K cvit )—2R^. (39) 



and the peak position occurs at: 

J ^ 22 - '0- 3- 2 )f-9-i (}-,}* 

4 

(72~7i)(6-3t2) 

\ 2(2 ~ 7 2) _(4-3 72 )(9-67 2 ) 

1 'l + 6 2 Ar) 2 . 



M peak = ^ c T^exp ^ uw 7 ua; " (40) 



In any case, the peak position should be close to the maximum of the two values given by 
eqns. ( 13T|) and ( HOI , respectively: 

(72~7i)(6-3-y 2 ) 

I ( P \ 2(2_72) /_ t2 , ,2-. ~( 4 -372)(9~672) ,„ .n . ,o, -(4-3 7 i)(9-6 7 i) . 

M peak ~ max ( ( (1 + 6 2 >1 2 ) 4 , (1 + 6 2 ^ 2 ) 4 | (41) 

Figure [5] shows the peak position as a function of exp(cx 2 ) = 1 + b 2 Ai 2 . The solid lines 
display the peak position (eqn. (j4"01 ) for 72 = 1.1 and for various values of p crtt /p =1, 10, 10 2 
and 10 3 . The dashed line corresponds to the peak position for 71 = 0.7 (eqn. f3"TI) whereas the 
dotted line is for 71 = 1 (isothermal gas). For realistic values of exp(cr 2 ) ~ 5 — 50, the peak 
position is almost always given by the solid lines since they lead to values larger than the 
values displayed by the dashed line. Therefore, the peak position is essentially determined 
by the nearly isothermal regime (effective polytropic exponent 72 ~ 1.1 in this case). This 
stems from the fact that as shown in fig. El the peak position changes very rapidly with the 
value of 7. 

The comparison between the dotted and solid lines reveals that, for realistic values 
of the ratio p crit /p ~ 10 2 — 10 3 , the peak position occurs at much smaller masses in the 
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Table 1: Parameters of clouds considered in fig. [HJ L, the size of the cloud and A , the pseudo 
Jeans length, are in pc, d is defined by eqn. P2|) . n, the average cloud density is in cm -3 , M d 
the cloud mass is in solar mass whereas V rmsi the root mean square velocity is in km s -1 . The 
parameters A4, (A4*) 2 and L/Xq which directly enter the theory, are dimensionless. Finally 
E k /Ep, E th /E p represent the corresponding ratios of kinetic and thermal over gravitational 
energy. 
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barotropic case considered here than in the isothermal case. This indicates that the cloud's 
mean density has to be much smaller in the barotropic case than in the isothermal case for 
the peak of the CMF to occur at the same mass. 

It is already interesting at this stage (detailed comparisons are performed in § 16. ip . to 
compare these analytical results with the numerical results obtained by Jappsen et al. (2005), 
displayed in their figure 5, which portrays the IMF obtained for various values of n crit /n ~ 0.5, 
5, 50 and 500, corresponding to n ~ 8.4 x 10 4 cm" 3 and n crit = 4.3 x 10 4 , 4.3 x 10 5 , 4.3 x 10 6 
and 4.3 x 10 7 cm -3 , respectively. We estimate the peak of the IMF for these four values 
of n cnt /n to be located at respectively M pea k/M & ~ 2.5, 0.5, 0.16 and 0.06, in solar mass 
units. These numbers can be compared with eqn. fj4*0l which, for 71 = 0.7 and 72 = 1.1, 
predicts that M peak oc (p C rit/p) ' 6 - Computing the corresponding exponent from the above 
quoted values, we find typically M pea k oc (p C rit/p) ' 55 in the simulations which, given the large 
uncertainties inherent to such difficult numerical calculations, appears to be fully compatible 
with the value of 0.6 predicted by eqn. (1301) . 

5.3. Results for a barotropic equation of state 

Figure M shows the core mass function for M. = 6 and hA\ = 2, for various values of 
the constant p crtt /p, namely 10, 10 2 and 10 3 . The solid lines display the case m — 1 in 
eqn. (l34"P . the dashed lines correspond to m = 3 whereas the dotted line corresponds to the 
isothermal case. As expected, the number of high-mass stars does not depend on the gas 
thermodynamics whereas the number of low-mass objects increases with increasing values 
of p crtt /p. For p cnt /p ~ 10, 71 = 0.7 and 72 = 1.1, the CMF is close to the isothermal case. 
Figure El also shows the non negligible dependence of the CMF upon the parameter m, which 
describes the - ill defined - transition between the two cooling regimes. 

6. Comparison with numerical simulations and observations 

In this section, we carry out detailed comparisons between our theory and the IMF 
determined both by the observations and by the numerical simulations performed by Jappsen 
et al. (2005) who determine the mass spectrum of collapsed objects in their simulations. 
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Fig. 7. — Comparison between our analytical IMF and the results of Jappsen et al. (2005) 
obtained using SPH simulations and sink particles, for the same characteristic parameters 
(see text). Top panel: n cnt = 4.3 x 10 6 cm" 3 ; bottom panel: n cnt = 4.3 x 10 7 cm" 3 . 
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Fig. 8. — Comparison between the Chabrier's IMF (dashed line) and the CMF predicted 
in the present theory for clouds following Larson-type relations (see text), for uq — 1 km/s 
(see eqn. ( 142ft ). The dash-curve which peaks at about 1M S is the Chabrier's IMF shifted 
upward by a factor 3 in mass to account for the shift between the CMF and the IMF. In the 
bottom-left panel the second dashed curve shows the original (non shifted) Chabrier's IMF. 
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Fig. 9. — Comparison between the exact CMF/IMF and the approximated one (given by 
eqns. ( ]A4j) -( lA3l) . The solid line and dashed lines show respectively the exact calculation for 
m = 1 and m = 3 while the dotted line shows the approximated CMF/IMF. 
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6.1. Comparison with numerical simulations 

In order to compare the simulations of Jappsen et al. (2005) with the present theory, 
we need to estimate the values of the parameters M., AA* and K CTlt in the simulations. Since 
the mean density in the simulations is 8.4 x 10 4 cm -3 while C s = 0.2 km _1 , the pseudo 
Jeans length (recalling that the pseudo Jeans length differs from the standard definition by 
a constant factor and assuming aj = 7r 5//2 /6) is about Aj = (y/n/2)C s / \JGp ~ 0.042 pc. 
The box size is 0.29 pc, giving L/\°j = 6.9. The Mach number in the simulation is 3.2 which 
gives .M* = -j=M.(\j/ L) v ~ 0.8. Finally, Jappsen et al. (2005) explore various values of 
n cnt . For n cnt = 4.3 x 10 6 cm -3 (fourth row of fig. 5 of Jappsen et al. 2005), this yields 
p crit / p ~ 50 while for n cnt = 4.3 x 10 7 cm -3 (fifth row of their fig. 5), this yields p CT \t/p — 500. 
Finally, we adopt the value m = 3 which leads to a temperature distribution very close to 
the one adopted by Jappsen et al. (2005) and a value of b = 0.5 typical of the type of forcing 
used in these simulations. We stress that our goal in this section is not to determine the 
most accurate IMF but rather to test the predictions of the analytical model with the direct 
numerical simulations. 

Dotted lines of fig. [7] show the numerical results of Jappsen et al. (2005) for the cases 
n cnt = 4.3 x 10 6 cm -3 (top panel) and n crtt = 4.3 x 10 7 cm -3 (bottom panel) while solid 
lines display the prediction of the present theory with the aforementioned parameters. The 
CMF predicted by the theory has been arbitrarily normalized to the numerical one near 1 
M Q . No adjustment has been made along the x-axis, implying that the peak position results 
from the dimensionless numbers (mainly Mach number) and the Jeans mass, including the 
coefficient aj. The agreement between the analytical theory and the numerical calculations 
is encouraging and suggests that the theory is indeed grasping the essential of the physics 
occurring in the simulations. In particular, it can be seen that in both cases, the peaks are 
approximately at the same place. It seems, however, that the analytical model is producing 
more brown dwarfs than the simulations. A possible explanation is that the value of aj 
should be modified. Indeed, aj = 7r 5 / 2 /6 is obtained by computing the mass enclosed within 
a sphere whose diameter is the Jeans mass. Although reasonable, such a definition is not very 
accurate and the comparisons suggest that we should take a value for aj about twice larger. 
Other, not exclusive, possibilities are that either the theory or the numerical simulations are 
not accurate enough in this regime. 

While the slopes at high masses are similar, it is maybe the case that the analytical 
curve never ceases to have a curved shape while the distribution inferred from the numerical 
simulations display a clearer powerlaw behaviour, particularly for the second panel. We note 
however that the numerical simulations appear to be a bit noisy, preventing to reach a definite 
conclusion at this stage and this behaviour has to be confirmed by further simulations. One 



-26- 



possibility would be that a more accurate density PDF should be used or that accretion onto 
the sink particles modify the shape of the distribution. 

6.2. Comparison with observations 

6.2.1. Cloud and model parameters 

Using the barotropic eos specified in § EJ our fiducial model is defined as, C s ~ 0.26 km 
s _1 for a density of n — 10 4 cm -3 (implying a temperature of about 20 K at this density), 
and m = 2 (see eqn. fl34|) ). n cnt = 4 x 10 5 cm -3 , ^ = 0.7, 72 = 1.07, and setting b = 0.75 and 
1] = 0.45. The theory entails three main parameters, p, A4, and Ai*, respectively the average 
cloud density, the Mach number at the cloud size and the Mach number at the mean Jeans 
length. Note that p enters the constant K cvlt stated by eqn. f|5B|) and also the normalization, 
in the definition of Mj. Finally, as mentioned previously, we choose aj = 7r 5 / 2 /6. 

Three parameters still represent a significant range of possibilities which are not all worth 
exploring, so we will restrict ourselves to a smaller subset. Since it is well established that in- 
terstellar clouds seen in the CO lines follow Larson-type relations, it seems natural to consider 
only clouds whose characteristics obey such relations. As recalled in § [2], the Larson relations 
(Larson 1981) correspond to V Tms ~ 0.8 kms -1 (L/lpc) 77 , fi~3x 10 3 cm _3 (L/lpc)~ 1 , where 
i] ~ 0.4 — 0.5. Therefore, once the size of the cloud is specified, all the relevant parameters 
can be inferred. Such a size-velocity dispersion scaling relationship has been confirmed by 
several authors (e.g. Heithausen et al. 1998, Falgarone et al. 2004) in various regions, 
although a value —0.7 in the size-density relation, instead of —1, seems to provide a better 
fit to the data. The mean density at a given size, however, appears to undergo significant 
variations. For example, the data compilation presented in Falgarone et al. (2004) seems to 
indicate that the typical density at a given size can vary from one region to another over at 
least one order of magnitude. This leads us to consider clouds of mean density defined as: 

n = (d x 10 3 cm" 3 ) ( — \ , V ims = (u x 0.8 km s" 1 ) ( — } . (42) 
Vlpc/ \lpc/ 

Note that theoretically, the existence of a density-size relation has been questioned by Kegel 
(1989), Vazquez- Semadeni et al. (1997) and Ballesteros-Paredes & MacLow (2002). More 
recent works have nevertheless obtained clear density-size relation (Hennebelle et al. 2007, 
Kritsuk et al. 2007). In any case, this relation is not at all fundamental in the present context 
and simply used as a simple prescription to determine the non dimensional parameters which 
control our theory. 

In the following, we will consider four different cloud sizes, namely 0.5, 1, 2 and 5 pc 
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representative of the size of a clump in a GMC We restrict ourselves to 3 values of do, namely 
3, 10 and 15, which correspond to clouds respectively as dense as and denser than the value 
originally quoted by Larson. Since it is interesting to explore also the dependence of our 
results upon the velocity dispersion, we also consider the case of various velocity dispersions, 
namely uq = 1, 1.5 and 2. Indeed, the recent study of Heyer et al. (2008) suggests that the 
velocity dispersion is correlated with the column density of the cloud E as w oc E 1 / 2 . Tabled] 
lists the various physical parameters of our cloud models, the dimensionless parameters which 
enter in the theory, as well as the ratios E^/E p and E^fE v where E t ^ = (3/2)MC 2 is the 
thermal energy, E k = (l/2)MV r 2 ms the kinetic energy and E p = (3/5)GM 2 /(L/2) is the 
gravitational energy. 



6.2.2. Model versus observations 

Figure M portrays the comparison between the models (solid lines) described in the 
previous section and the Chabrier's system IMF (dashed lines, Chabrier 2005). Indeed, we 
recall that the present theory predicts the distribution of prestellar cores, as observed in 
dust-continuum surveys, and does not include further subfragmentation of these cores into 
individual objects (see discussion in §7.1.2 of paper I). In the same vein, the resolution of 
present observational surveys does not allow to resolve individual objects in the observed 
bound cores. Since there is observational evidence that the CMF appears to be shifted with 
respect to the IMF by a constant factor of about 3 in mass (Motte et al. 1998, Alves et al. 
2007), the Chabrier's IMF is shifted in mass by the same factor on the figures. The Chabrier's 
IMF is normalized to the observed stellar density near 1 M in the solar neighborhood 
but the normalization in our calculations is arbitrary. Therefore, in order to compare the 
observational and theoretical curves, they have been adjusted along the y-axis at the same 
peak value. We emphasize that no adjustment has been made along the x-axis (mass). That 
is, for a given analytic distribution J\f(A4), we plot dN/dlogM = ln(l0)MAf(Ai), where 
M = M°jM. 

For do = 3, which corresponds to the value originally quoted by Larson, the agreement 
between the theoretical CMF and the Chabrier's IMF is poor. The two main important 
disagreements are the peak value, which is larger for the model than for the shifted Chabrier's 
IMF, and the slope at large masses which is not stiff enough for the model. While the first 
discrepancy is simply due to the Jeans mass being too large, the second discrepancy can be 
understood as follows. Equation (40) of paper I shows that at large masses, the value of x 
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is slightly more complex than indicated by eqn. (42) of paper I, and is exactly given by: 



Although the value of the second term is small when .M* — 1 and Ai ~ 6 — 10, it can lead 
to non negligible corrections when M.* takes larger values and A4 smaller values, as is the 
case here for some of the cloud conditions listed in table 1, for which the second term can 
amount to about 0.2, leading to x ~ 1.1. 

For do = 10 and uq = 1, a better although imperfect agreement between the observations 
and the prediction of the present theory is obtained. We see that the peak position, even 
though closer to the peak of the shifted Chabrier's IMF, does not exactly coincide with 
it. The slope at large masses is very close to the Salpeter value, x ~ 1.35, but too many 
massive stars are produced. There is also an obvious disagreement in the low-mass regime 
where the model predicts significantly less brown dwarfs and low-mass stars. The agreement 
is improved if we consider a higher velocity dispersion, uq = 1.5, as revealed by the first 
column of fig. [101 The peak is now approximately at the right place whereas the number 
of high mass stars is in better agreement, although still slightly larger than the observed 
values. As noted in §7.1.4 of paper I, the peak of the IMF does not change significantly 
when changing the size of the clumps, assuming they obey Larson-type relations. This is 
due to a cancellation between the increasing Jeans mass and the increasing Mach number as 
the cloud size increases. 

In the case do = 15, a better agreement between the theory and the shifted Chabrier's 
IMF is obtained both for uq — 1 (fig. M) and uo = 1.5 (fig. [TO]) . The peak and the high 
mass part of the distribution are very well reproduced by the models except perhaps for 
L = 0.5 pc and Uo = 1, for which the number of high mass stars may drop too stiffly 
(because of the too low global Mach number A4). The lack of brown dwarfs is less severe 
than for the case do = 10. For the cases L = 2 and 5 pc and uq = 1.5, in particular, the 
remaining disagreement is occurring at very low masses where the observations are not well 
constrained. Finally, in order to investigate the influence of larger values of Uq, we also 
present (third column of fig. [TU]) the case do = 15 and Uo = 2. In that case, the peak occurs 
at too small masses while the slope in the high mass part is slightly too shallow. 



x = 



n + 1 6(1 - r])\nM* 
2n-A ~ 277 + I 



(43) 



6.3. 



Discussion 



According to the predictions of our theory, we conclude that typical conditions for 
star formation should differ from those specified by the Larson's relations. In particular, 
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the clouds characteristic of star forming regions need to be denser than predicted by these 
relations. However, two points must be stressed. First of all, while comparing our results 
with the shifted Chabrier's IMF (based on the observed shift between the CMF and the 
IMF), we choose a constant value of 3 for the shifting factor. Although such a value has 
been inferred from observations in the Pipe nebula (Alfves et al. 2007), it may in fact 
depend on the cloud and exhibit some variation. For instance, Motte & Andre (1998) find 
a value closer to 2 in the rho Ophiuci cloud, while J0rgensen et al. (2007) find a ~ 10-15% 
core-to-star efficiency in Perseus, implying a ~ 7 to 10 shifting factor between the CMF 
and the IMF, although such numbers should be taken with due caution. Second of all, our 
theory ignores the magnetic support, which would add up to thermal and turbulent support. 
As discussed in §4.4 of paper I, magnetic support can enter our theory by simply rescaling 
the turbulent support, and thus the effective Jeans length/mass, according to the observed 
correlation between the magnetic intensity and the density and rms velocity dispersion in 
molecular clouds, affecting also to some extend the Mach number and thus the width of the 
lognormal distribution. However, such a treatment is certainly a crude approximation of 
the true impact of magnetic field. At any rate, taking magnetic field into account would 
increase the support, implying denser clumps. All these points suggest that the clumps in 
which star formation is taking place should be at least as dense as, probably denser than the 
do = 15 case displayed in the previous section. This suggests that star formation should occur 
preferentially in dense gas regions where the gas density is typically 3-5 times larger than the 
canonical values originally inferred by Larson for the CO clumps. Note that this must be the 
case at all scales which, in principle, arises in a turbulent flow. Indeed, polytropic flows are 
scale free and varying the mean density implies varying the density at all scales. For a clump 
of 1 pc, the density for which we get reasonable agreement is of the order of > 1.5 x 10 4 cm -3 
according to our analysis. This may suggest that a triggering mechanism might be necessary 
to lead to such high mean densities. This is consistent with a picture of star formation 
happening in relatively dense clumps, as suggested by many observations of star forming 
regions. In Perseus, Serpens, Ophiucus, for instance, cores are found preferentially at high 
cloud column densities (Av > 6), corresponding to mean densities n > 2 x 10 4 cm -3 (Enoch 
et al. 2008). These conclusions are also consistent with the low star formation efficiency 
observed in the Galaxy since they imply that only a tiny fraction of the gas is actually 
forming stars efficiently. 

There is a possible trend for the models to form too few low-mass brown dwarfs. Al- 
though quantifying the exact amount of this discrepancy requires more sound determinations 
of the IMF in the low-mass brown dwarf regime, its origin, if true, is unclear. One possibility 
is that the thermodynamics of the gas is different from what has been assumed here. In par- 
ticular, if the critical density, p crit , or the effective adiabatic exponent, 72, are respectively 
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higher and/or smaller than the value used in the present study, the number of small mass 
objects will increase. Since the thermodynamic properties of the high density gas are not 
well known this possibility can not be excluded. Another suggestion is that gravity, which 
appears in our theory only under the form of a Jeans criterion to define gravitationally 
bound objects, might lead to more fragmentation in the very-low-mass regime. Similarly, a 
lognormal distribution for the turbulence-driven density fluctuations and a variance given 
by eqn.fjS]), in particular in the non-isothermal and/or magnetized case, might be an over- 
simplification. As mentioned earlier, simulations including compressive forcing, for instance, 
clearly show a departure from such a distribution, with a PDF extending further down into 
the small-scale regime (Federrath et al. 2008). In the same way, the density PDF that we 
consider is valid in the isothermal case but should clearly be modified when non isothermal 
eos are considered. Let us recall that Passot & Vazquez-Semadeni (1998), find that the PDF 
is not a lognormal and present a powerlaw tail at high densities when 7 < 1. Taking this 
effect into account would certainly increase the number of brown dwarfs. Finally, it sounds 
likely that gravity will also lead to a similar effect with the high density tail of the PDF 
becoming less stiff and more extended (Klessen 2000, Vazquez-Semadeni et al. 2008, Hen- 
nebelle et al. 2008). Indeed, the role of gravitational focussing in the production of brown 
dwarfs has been explicitly investigated by Bonnell et al. (2008). 

The global success of the present theory reinforces the suggestion that gravoturbulent 
fragmentation seems to be the dominant global physical mechanism responsible for star and 
brown dwarf formation. This is supported, on the theoretical side, by the general fairly good 
quantitative agreement, over basically the entire mass spectrum, between the observationally 
determined and the theoretical CMF/IMF, as illustrated in the previous sections, and on the 
observational side by various explorations of young brown dwarf properties (Luhman et al. 
2007, Duchene et al. 2009) and by the recent statistical determination of the brown-dwarfs 
to stars ratio in several young clusters (Andersen et al. 2008). This latter study shows 
that a Chabrier IMF, continuously extending from the stellar to the brown dwarf regime, is 
fully consistent with the observations. It is not excluded, however, that the low mass end 
of the IMF might be affected by secondary mechanisms such as, for instance, dynamical 
ejections and radiative feedback (Bate 2009) or massive disk fragmentation (Stamatellos & 
Whitworth 2009). Concerning the second mechanism, it should be stressed, however, that 
these calculations do not include magnetic field, which has been shown to severely inhibit 
disk fragmentation (Machida et al. 2005, 2008, Hennebelle & Teyssier 2008, Price & Bate 
2008) or even to prevent disk formation by magnetic braking (Galli et al. 2006, Price & 
Bate 2007, Hennebelle & Fromang 2008, Mellon & Li 2008a,b). Interestingly, Machida et 
al. (2008) find that the formation of low mass objects is possible during the second collapse 
after the gas and the magnetic field decouple. 



7. Conclusion 



In this paper, we have further investigated the analytical theory we have recently devel- 
oped for turbulence-driven star formation (Hennebelle & Chabrier 2008). We have explored 
and quantified the influence of various properties of the flow, such as the turbulent forcing or 
the fluctuations of the rms velocity field, on the mass spectrum. Moreover, we have extended 
our previous calculations, based on an isothermal equation of state for the gas, to a more 
general polytropic or barotropic case. We show that the low mass end of the CMF/IMF 
is very sensitive to the value of 7, the effective adiabatic index, and to the temperature- 
density dependence within the cloud (see Fig. 4 and 6). Departure from isothermality with 
smaller values of 7 (7 < 1) leads to a more compressible flow and cooler overdense structures, 
increasing the number of small-mass collapsing cores, precursors of brown dwarfs. Exten- 
sive comparisons with the non-isothermal numerical simulations performed by Jappsen et 
al. (2005) show the good agreement between the IMF obtained in the simulations and the 
predictions of the theory. 

Taking into account the thermal properties of the flow and a more complete description 
of turbulence forcing improves the quantitative agreement between the theoretical CMF/IMF 
and the observationally-derived IMF, for more realistic Mach number values, typical of star 
forming regions. A proper combination of these two properties, through the 7 and b coeffi- 
cients, leads to an excellent agreement between theory and observation. Unfortunately, the 
lack of knowledge of both the precise thermodynamic properties of the flow at high density 
and of the PDF of non-isothermal turbulence prevent an accurate determination of these two 
parameters under star-forming cloud conditions. It is not excluded, in particular, that the 
theory shows a trend for a deficiency of brown dwarfs. As mentioned above and in paper I, 
however, observational determinations of the brown dwarf IMF in the low-mass regime are 
still very uncertain and more accurate determinations are required before allowing robust 
conclusions. Various comparisons between our models and the observed IMF, for different 
cloud conditions, suggest that star formation should dominantly occur in clouds typically 
5 times denser than the CO clumps characterized by Larson (1981), a point which seems 
to be supported by observations. This raises the question of the universality of the IMF 
since choosing different cloud parameters would lead to different CMF/IMF, in contradic- 
tion with the observational evidence. The most plausible answer to this apparent puzzle 
lies in the universal nature of the turbulence spectrum, which tends to produce clouds with 
similar characteristics (mass versus size relation, velocity dispersion...), with efficient star 
formation, i.e. formation of gravitationally bound cores, only or preferentially occurring in 
the clouds exceeding some mean density value. Note that in the present theory, the char- 
acteristic mass of the IMF/CMF is set up by the cloud's general conditions. Meaning it 
depends both on the average density and temperature (thus mean thermal Jeans mass) of 
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the cloud but also on its characteristic (large-scale) Mach number, illustrating the combined 
effect of gravity and turbulence to set up the IMF. However, as demonstrated in §7.1.4 of 
paper I and discussed in §5.2 of the present paper, the similar but opposite dependence of 
the Jeans mass and Mach number upon the cloud's size/mass, for clouds obeying Larson's 
type relations, leads to an IMF peak position and characteristic mass, which arise in the 
nearly isothermal regime (see §5.2), nearly independent of the cloud's conditions. This is in 
contrast to the suggestions of Elmegreen et al. (2008) and Bate (2009). While the former 
authors invoke the weak environmental dependence of the ocurence of gas-dust coupling to 
explain the universality of the characteristic Jeans mass, the second author argues for radia- 
tive feedback of prestellar cores to yield an effective Jeans mass nearly independent of the 
cloud's density. Interestingly, these two approaches yield opposite density-dependences of 
the Jeans mass. More detailed comparisons between dedicated simulations and the present 
theory are necessary to explore this "universality" issue. 

An already interesting success of the present theory, testable with numerical simulations, 
is that it provides a consistent explanation to why gravitational clumps obtained by taking 
into account purely thermal support lead to too steep an IMF at large masses (Padoan et 
al. 2007), while simulations which use sink particles and include by definition the turbulent 
support with a Kolmogorov-like power spectrum index obtain the proper Salpeter slope (e.g. 
Klessen 2001, Bate & Bonnell 2005). Furthermore, as explored in §6.3 of paper I, this theory 
shows that turbulence globally decreases star formation efficiency, explaining at least partly 
the low efficiency of star formation observed in the Galaxy. 

The overall good quantitative agreement between the present results and the CMF/IMF 
derived both from numerical simulations (including compressible-mode forcing and non- 
isothermal gas eos) and observations over basically the entire mass spectrum, down to about 
0.01 M Q , although by no means an irrefutable proof, reinforces the validity of the present 
theory. This latter suggests that gravo-turbulent fragmentation, including the two opposite 
roles of turbulence, turbulence-driven density enhancements at the cloud's scale on one 
hand and turbulent support at the core's scale on the other hand, is a promising dominant 
physical mechanism responsible for star and brown dwarf formation. Alternative mechanisms 
very likely contribute only marginally to this dominant process at the stage of prestellar 
core formation. For instance, observations (Andre et al. 2007, Evans et al. 2008) suggest 
that collision timescales between cores are substantially longer than their lifetimes, so that 
collisions and thus competitive accretion between dense cores are unlikely to be dominant. 
This in turn suggests that the prestellar dense cores are likely to evolve individually into 
stellar systems. This property lies at the very heart of the present theory. As mentioned 
earlier, various observations of star forming regions (Motte & Andre 1998, Andre et al. 2007, 
Alfves et al. 2007, Enoch et al. 2008) confirm the similarity of the CMF and IMF, providing 
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further evidence that stellar masses are directly linked to the core formation process. In this 
picture, at the heart of the present theory, the total (including further multiple systems) 
stellar masses are determined by the initial fragmentation of the cloud/clump into dense 
cores. This is in contrast to alternative mechanisms where the stellar masses - thus the 
IMF - are dominantly determined by gas-to-star conversion processes, such as competitive 
accretion, feedbacks or winds, for which the final IMF is not expected to reflect the CMF. 

Understanding the transformation of the prestellar core mass function (CMF), as de- 
scribed by the present theory, into the stellar (system) mass function (IMF) remains an 
open issue. The aforementioned observations suggest a nearly mass-independent conversion 
efficiency factor, in the ~ 10-50% range, at least within the sensitivities of present sur- 
veys. Such a nearly constant star formation efficiency of cores is predicted theoretically if 
magnetically-driven outflows are the primary mediating factor from cores to stars (Matzner 
& McKee 2000). Numerical experiments, on the other hand, seem to suggest that the shape 
of the CMF is barely affected by the nature of the core-to-star, i.e. CMF to IMF, conversion, 
except possibly in the very brown dwarf regime (Swift & Williams 2008). Again, a more 
precise determination of the IMF in this domain is needed to further explore this issue. 

A substantial fraction of the prestellar cores, however, will fragment into multiple in- 
dividual objects. This sub-fragmentation episode very likely involves strong gravitational 
interactions and angular momentum redistribution between the individual objects within 
a dense core. Processes like radiative feedback, dynamical ejections, disk fragmentation, 
magnetic outflows might play some important role at this stage. An analytical description 
of these strongly non-linear mechanisms seems to be rather elusive and we emphasize that 
the present theory only addresses the initial stages of star formation, i.e. the formation of 
gravitationally bound prestellar cores. As such, this theory provides a sound theoretical 
foundation for the formation of pre- stellar /brown dwarf cores from the turbulent fragmenta- 
tion of a cloud, and the related CMF/IMF. We also stress that, although relying on the same 
general concept of gravo-turbulent fragmentation as the Padoan & Nordlund (2002) theory, 
the present theory differs both qualitatively and quantitatively from the Padoan-Nordlund 
one, with testable predictions, as discussed in §7.2 of paper I. Although the exact mass dis- 
tribution derived in the present theory in the general (non-isothermal) case requires a (very 
simple) numerical resolution (eqn. fl26l) ). we provide in Appendix A a reasonably accurate 
analytical approximation of this CMF/IMF. This analytical expression should be used as a 
benchmark to compare core/system mass distributions obtained from numerical simulations 
or inferred from observations, with the present theory. The HERSCHEL satellite, launched 
in spring 2009, should reveal very- low- mass gravitationally bound prestellar cores and probe 
the CMF well into the substellar domain in various star forming regions. This will provide 
a stringent test for the mass distributions and the prestellar cloud conditions predicted by 
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A. A simple and accurate approximation for the CMF/IMF 

As discussed in § 14. 2\ in the non isothermal case the criterion stated by eqn. (1261) 
cannot be obtained analytically. Since this makes the calculation of the CMF/IMF less 
straightforward, we propose in this appendix an accurate and simple approximation. 

In the purely thermal case, the criterion for collapse can easily be obtained for any value 
of 7 and is given by eqns. (1241 and fl26l) . for V rms = and .M* = 0. Our approximation 
consists to approximate the condition (T2"4"l) by: 

M > Mf + Mf rb , (Al) 

i.e. to assume that in order for the fluid cell to be gravitationally unstable, its mass must 
be greater than the sum of the critical masses corresponding to each individual support, 
thermal and turbulent. In the case of a barotropic equation of state, as the one considered 
in § [5j this gives for the critical mass: 



3 71 -4 37 2 -4 2 +1 

M > Mg = M ?(71 ) (^) + M S '( 72 ) (^) + M°' M Ml (^) (A2) 



This leads to: 



_ 3 71 -4 1 _ 37 2 -4 

M > M C R ~ K 7 ^ + K^ 2 + MlR 2v+1 . (A3) 



where R = R/X°j(^i), M = M/M°'( 7l ) and K ciit = (p crit /p)^-^. Thus 

= ?iLzi^-. + K ^j2^im-i + Ml[2n + (A4) 

ctit 71 — 2 72 — 2 
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These quantities enter the final CMF/IMF, given by eqn. fl27|) . that we rewrite here: 
»r/^ P 1 dR dd c R 1 / (<^) 2 5% a 2 \ 



(A6) 



6% 






_dM R 


dR 


dR 



We also recall that 



i v f\°; 



a 2 ~ln(l + &lM 2 ), (A7) 

where Aj, given by eqn. (122|) . reflects the velocity dispersion at the Jeans length characteristic 
of the cloud's average density. We now have an explicit set of equations straightforward to 
implement. A comparison between this approximate expression and the exact one (solving 
eqns. (jUj) and fl26l) ) for the second case displayed in fig. [6], is shown in fig. [9j As can be seen, 
the agreement is fairly reasonable, in particular the case m = 1 (see eqn fl34|) for definition). 



This approximated CMF/IMF should be sufficient for a first quick-and-easy comparison 
between the present theory and either observational or numerical results. More accurate 
comparisons, however, necessitate the proper resolution of eqn. (1261) . a fairly easy task. 



B. General formulation of the CMF valid for any equation of state 

In this appendix, we write the equations to be solved in the case of a general equation 
of state written as: 

~ f(p) = /(% (Bl) 



Proceeding as before, we obtain: 



M c R = R[f(j^)+MlR 2 r>), (132) 



-36- 



leading to 



dM% /(§)-3§f(§) + (2r7 +lVW^ 



dR i Lp 

which generalize eqns. (1A3|) and (1A4p while eqns. (1A5[) - (IA7[) remain unchanged. As described 
in § eqn. (lB2p must be solved numerically. 



C. Core lifetime and the IMF 

The relationship between the core mass function and the initial mass function has been 
questioned by Clark et al. (2007) on the basis of a timescale argument. 

Their point is the following: since the cores are dynamically collapsing, their typical 
timescale of evolution is the freefall time, r// ^ O n the other hand, it is expected 

that the mass of the core should be about one Jeans mass, M ~ Mj oc C^/yfp. Thus, 
M oc Tff. This implies that while a core of mass M is collapsing, about 10 cores of mass 
M/10 should have time to form and collapse. Therefore, in order to produce an IMF having 
the slope of Salpeter, x = 1.35, the CMF should have a slope x = 1.35 — 1 = 0.35, not 
compatible with the observations. 

We suggest that the solution to this problem relies on the fact that what determines 
the mass of the massive cores is not the thermal Jeans mass but the turbulent Jeans mass. 
Indeed, as discussed in paper I, this is the turbulent dispersion (taken into account in the 
parameter M.*) which, in our theory, leads to the slope of Salpeter, x = 1.35. For the 
turbulent Jeans mass, we have Mj urb oc Vf ms j y/p. 

Thus, we have: 



R 3v 

Mcx^-, (CI) 
\ff> 



and since M oc pR 



-l + i? 1-17 

p oc MW,^ oc M^. (C2) 

Taking, for example r\ = 0.5, we have Tff oc M 1 / 4 to be compared with the relation 
Tff oc M obtained by Clark et al. (2007). If, as it is the case in the present theory, the mass 
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of the massive cores is determined by the turbulent dispersion, their density varies with their 
mass in a much shallower way than if their mass would solely be determined by the thermal 
support. 

The implication is that the CMF has to be only slightly shallower than the IMF in 
order to obtain the Salpeter's slope. Recalling that the slope given by eq. (H3l) is indeed 
slightly shallower than the Salpeter value, we believe that this is fully consistent with the 
core mass function being at the origin of the IMF. A more quantitative estimate including a 
full calculation is beyond the scope of this paper and will be considered elsewhere. 
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Fig. 10. — Same as fig. [8]for do = 10 and 15 and for various values of L and uq (=1.5 and 2). 



